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^. ■ ABSTRACT 

Q^l Wc have performed a Monte-Calro simulation for Galactic population of pul- 

O , sars and for the 7-ray observations. Wc apply two-layer outer gap model, 

H ■ 

i which has been developed by Wang, Takata & Cheng, for the 7-ray emission 

process, and study the radiation characteristics as a function of the magnetic 
inclination angle (a) and the Earth viewing angle (C). In our model, the 7-ray 
flux and the spectral cut-off energy tend to decrease as the inclination and 
viewing angles deviate from 90°. The emerging spectrum above 100 MeV be- 
■ comes soft with a photon index p~^1.8 — 2for<;^— > 90° and p ^ 1.2 — 1.3 



for Q <C 90°. Our simulation predicts that the pulsars with larger inclination 
angles (a = 70 — 90°) and larger viewing angles = 70 — 90°) have been 
preferentially detected by the Fermi 7-ray telescope, and hence the observed 
pulse profiles of the 7-ray pulsars have the double peak structure rather than 
single peak. In the simulation, most 7-ray millisecond pulsars are categorized 
as the radio-quiet 7-ray pulsars, because its radio fluxes are under the sen- 
sitivities of the major radio surveys. Even we drastically increase the radio 
sensitivity by a factor of ten, the number of radio-selected millisecond pul- 
sars detected by the Fermi ten years observations is still much less than the 
expected 7-ray-selected millisecond pulsars, indicating the radio-quiet mil- 
lisecond pulsars must contribute to the Fermi unidentified sources and/or 
the 7-ray background radiations. We argue that 7-ray pulsars observed with 
a smaller viewing angle (C ^ 90°) will appear as low-efficient 7-ray pulsars. 
For example, unique radiation properties of the low-efficient 7-ray pulsar, 
PSR J0659+1414, can be explained by the present our gap model with a 
viewing geometry oi a ^ C, ~ 40° — 50°. 
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1 INTRODUCTION 

The Large Area Telescope [LAT) on board the Fermi 7-ray telescope has increased number 
of 7-ray pulsars, and the recent Fermi catalog includes more that 60 7-ray pulsars, which 

Saz Parkinson et al. 2010). 

Furthermore, the detection of radio millisecond pulsars associated with about 20 unidentified 
Fermi point sources (e.g. Ray 2010; Caraveo 2010; Ranson et al. 2011; Keith et al. 2011) 
has been reported, suggesting that the millisecond pulsar, as well as the canonical pulsar, 
is one of the major Galactic 7-ray source. It can be expected that more 7-ray pulsars will 
be added to the list over the Fermi mission. The spectral shape and the pulse morphology 
measured by the Fermi have been used to discriminate particle acceleration and 7-ray 
emission models; the polar cap model (Ruderman & Sutherland 1975; Daugherty & Harding 
1982, 1996), the slot gap model (Arons 1983; Muslimov & Harding 2004; Harding et al. 
2008; Harding & Muslimov 2011) and the outer gap model (Cheng, Ho & Ruderman 1986a,b; 
Hirotani 2008; Takata, Wang & Cheng 2010a). The polar cap model assumes the acceleration 
region near the stellar surface, and the slot gap and outer gap models assume the emission 
region extending to outer magnetosphere. The cut-off features of the 7-ray spectra of the 
Crab and the Vela pulsars measured by Fermi imply that the 7-ray emission site of the 
canonical pulsars is located in the outer magnetosphere rather than near polar cap region, 
which produces a cut-off feature steeper than the observed one (Aliu et al. 2009; Abdo et 
al. 2009c, 2010d). Romani & Watters (2010) and Watters & Romani (2011) have studied 
morphology of the pulse profiles of the young pulsars predicted by the our gap and slot gap 
models, and they argued statistically that the outer gap geometry is more consistent with 
the Fermi observations than the slot gap model. Venter, Harding & Guillemot (2009) found 
that the observed pulse profiles of several millisecond pulsars detected by the Fermi cannot 
be explained by the outer gap and/or the slot gap models, and proposed a pair-starved polar 
cap model, in which the particles are continuously accelerated up to high altitude because 
of the insufficient multiplicity of the pairs. 

The high quality data measured by the Fermi enables us to perform a detail study for 
population of the 7-ray pulsars. Takata et al. (2010a) have studied the relation between 
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the emission properties (luminosity and spectral cut-off energy) and pulsar characteristics 
(e.g. rotation period and magnetic field). They proposed that the outer gap accelerator 
model controlled by the magnetic pair-creation process can explain the observed population 
statistics better than that controlled by the photon-photon pair-creation process (Zhang & 
Cheng 1997, 2003). Wang, Takata & Cheng (2010) have fitted the observed phase-averaged 
spectra by using a two-layer outer gap model, in which the accelerator consists of a wide 
but low charge density main region and a narrow but high charge density screening region. 
They suggested that the relation between the 7-ray luminosity (L^) and the spin down 
power (Lsd) can be expressed as oc Lf^ with /3 ~ for Lgd ^ 10^^ erg/s, while f3 ~ 0.5 
for ^ 10^^ erg/s. This relation is consistent with the theoretical expectation, i.e. the 
gap fractional size / is determined by / = min{fzc,fm), where f^c and are the gap 
size determined by photon-photon pair-creation process (Zhang & Cheng 1997, 2003) and 
magnetic pair-creation process (Takata et al 2010a), respectively (c.f. section [3l2|) . Equating 
fzc and fm corresponds to Lsd ~ 10"^^ erg/s. 

The study for population synthesis of the radio pulsars has been developed by using the 
detailed modeling of the radio emissions and the radio surveys (e.g. Bailes & Kniffen 1992; 
Sturner and Demer 1996; Faucher-Giguere and Kaspi, 2006). Also the population study of the 
7-ray pulsars have been developed by several authors (e.g. Gonthier et al. (2002) for the polar 
cap model, Cheng & Zhang (1998) for the outer gap model, and Story, Gonthier & Harding 
(2007) for the slot gap model). For example. Story et al. (2007) studied the population 
of 7-ray millisecond pulsars with the slot gap accelerator model, and predicted the Fermi 
observations. They predicted that the Fermi will detect 12 radio- loud and 33-40 radio-quiet 
7-ray millisecond pulsars. With the Monte-Calro simulation of the outer gap, Takata, Wang 
& Cheng (2011a,b) have explained the observed distributions of the characteristics of the 
7-ray pulsars detected by the Fermi with the six-month long observations. They predicted 
that at least 80 7-ray canonical pulsars have a 7-ray flux exceeding the sensitivity of the 
six-month Fermi observations, suggesting the present observations have missed many 7- 
ray emitting pulsars. Watters & Romani (2011) simulated the Galactic distribution of the 
young 7-ray pulsars, and compared the simulated pulse morphology with the Fermi results. 
They argued statistically that the outer gap model explains the distributions of the pulse 
morphology (e.g. the phase-separation of the two peaks) measured by the Fermi better 
than the slot gap model. The population studies (e.g. Kaaret & Philip 1996; Cottam, Jean 
Faucher-Giguere & Loeb 2010; Takata et al. 2011b) have also pointed out that unidentified 
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pulsars, in particular the millisecond pulsars located at high-Galactic latitudes will associate 
with the Fermi unidentified sources (Abdo et al. 2010b), and will contribute to the 7-ray 
background radiations. 

In the our previous studies (Takata et al. 2011a,b), we ignored the dependence of the 
radiation characteristics on the viewing geometry (i.e. the magnetic inclination angle and 
the Earth viewing angle measured from the pulsar's rotation axis), and focused on the 
distributions of the properties (e.g. the rotation period and magnetic field) of the 7-ray 
pulsars. However, the observed radiation characteristics, such as the fiux, spectral cut-off 
energy and pulse profile, must be affected by the viewing geometry. To perform a more 
solid study on the population, which is compared with the high-quality Fermi data, it is 
required a three-dimensional model that takes into account the dependence of the radiation 
characteristics on the viewing geometry. 

In this paper, we develop a Monte-Carlo study for the population of the 7-ray pulsars 
with the radiation model that takes into account the dependence on the viewing geometry. 
In section [2], we review the Monte-Carlo simulation for the population of the 7-ray pulsars. 
We discuss our outer gap model in section [3l In section 14.11 we discuss the dependence of 
the radiation characteristics on the viewing geometry. In section 14.21 we compare the results 
of the Monte-Calro simulation with the Fermi six-month long observations. We also show 
the expected population of the 7-ray pulsars if the Fermi observations continue five-years 
or ten- years. In section [5l after summarizing our simulation results, we discuss the viewing 
geometry of PSR J0659+1414, which is known as a low efficient 7-ray pulsar. 



2 MONTE-CARLO SIMULATION FOR THE PULSAR POPULATION 

In this paper, we denote the canonical pulsar and the millisecond pulsar as CP and MSP, 
respectively. We assume that the birth rates of the CPs and the MSPs are ~ lO^^yr^^ 
and 10"^ ~ 10"^ yr~^ (Lorimer et al. 1995, Lorimer 2008), respectively, and we ignore 
the millisecond pulsars in globular clusters. The birth location is determined by the spatial 
distributions given by (Paczynski 1990), 

pr{R) = — ^2 ' 

Pz{Z) = -Le-l^l/^o.p^ (1) 

•^exp 
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where R is the axial distance from the axis through the Galactic centre perpendicular to 
the Galactic disk and Z is the distance from the Galactic disk, i?cxp = 4.5 kpc and an = 
[l-e-^"'''-/^'^-''{l + Rmax/Rexp)]~'^ with Rynax = 20 kpc. In addition, we apply Z^xp = 75 pc 
for the CPs and Z^xp = 200 pc for the MSPs, respectively. 

To obtain current position of each simulated pulsar, we solve the equation of motion 
from its birth to the current time. The equation of motion is given by 



dt^ R dR ' ^ ^ 



dZ^ ^ 

dt^ dZ ' ^ ^ 

and 

Rv^ = constant. (4) 

Here is the azimuthal component of the velocity, ^t^t = ^sph + ^dis + is the total 
gravitational potential, where ^sph, ^dis and $/i are spheroidal, disk and halo components 
of the Galactic gravitational potential, and are given by 

MR, z) = - , ™' (5) 
ir^ + K + {z^ + bfyi^? 



where i = sph and dis, asph = 0, bsph = 0.277 kpc, Mgph = 1.12 x IO^'^^Mq, adis = 3.7 kpc, 
bdis = 0.20 kpc, and M^is = 8.07 x 1O^°M0, while for the halo component 



r) = 



1 J 1^]^ _|_ '"^^ _|_ ^ -if''" 

2 \ rl I r \rr 



(6) 



where rc = 6.0 kpc and Mc = 5.0 x lO^^M© (c.f. Burton & Gordon 1978; Binney & Tremaine 
1987; Paczynski 1990). The Lagrangian in units of energy per unit mass is given by 

,.fMhI^.^^,(n,z), (7) 

where v is the velocity. 

For the initial velocity of each modeled pulsar, we assume random isotropic direction of 
the velocity, and assume magnitude drawn from a Maxwellian distribution with a charac- 
teristic width of ay = 265 km/s for the CP and ay = 70 km/s for the MSP (c.f. Hobbs et 
al. 2005), namely. 
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2.1 Pulsar characteristics 

In this section, we describe how we calculate the current value of the various characteristics 
of the simulated pulsars. 

2.1.1 Canonical pulsars 

For canonical 7-ray pulsar, which is born soon after the supernova explosion, the true age is 
more or less equal to the spin down age r = P/2P, where P and P are the rotation period 
and its time derivative. This allows us to calculate the present distributions of the pulsar 
characteristics from the initial distributions. 

For the Crab pulsar, the initial period is estimated to be Pq ~ 19ms. A short birth 
rotation period is also expected for the young 16ms pulsar PSR J0537-0691 (Marchall et al. 
1998), suggesting most pulsars were born with Pq ~ 20 — 30 ms. There is a good evidence 
that some pulsars were born with a longer initial period (e.g. Pq ~62 ms of PSR J191 1-1925, 
Kaspi et al. 2001). However, we found that the distribution of initial rotation period does 
not affect much to the results of the following Monte Carlo simulation. For example, we 
compared the results of two simulations with different distribution of the initial period; one 
is that all simulated pulsars are randomly distributed in the range Pq = 20 — 30 ms, and other 
is that 70 % of pulsars are distributed at Pq = 20 — 30 ms and 30 % are Pq = 20 — 100 ms. In 
such a case, we could not see any significant difference in the populations of the simulated 
7-ray pulsars. The difference in the simulated distributions of the 7-ray pulsars becomes 
significant if we assume that > 50 % of new born pulsars is distributed with the initial 
period of Pq ^ 30 ms. Because we do not know well the exact distribution of the initial 
period, we randomly choose the initial period in the narrow range Pq = 20 — 30 ms, which 
provides a more consistent result of our simulation with the Fermi observations. 

We assume a Gaussian distribution in log^gP^ for the initial distribution of the stellar 
magnetic field measured at the magnetic equator. 



In this study, we apply log^^o-^o = 12.6 and as = 0.1. Because the canonical 7-ray pulsars 
are younger than 10 Myr, we ignore evolution of the stellar magnetic field, which may be 
important for the neutron star with an age older than 10 Myr (Goldreich & Reisenegger 
1992; Hoyos, Reisenegger & Valdivia, 2008). 




(9) 
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With a constant stellar magnetic field in time, the period evolves as 

m=(n= + =^*) . (10) 

where Rg = 10^ cm is the stellar surface and / = lO^^gcm^ is the neutron star momentum 
of inertial. In addition, we artificially assumed that the magnetic inclination angle does not 
evolve with the spin down age, while a decreases on the spin down time scale has been 
pointed out (e.g. Davis & Goldstein 1970; Michel 1991; Tauris & Manchester 1998). The 
time derivative of the rotation period is calculated from 

2.1.2 Millisecond pulsar 

We assume that all MSPs are born through the so called recycled process, in which the 
accretion of the matter from the low mass companion star spins up the neutron star. It 
implies that the true age of the binary system is different from the spin down age of the MSP. 
However, we expect that the Galactic distribution does not depend on the spin down age of 
the MSPs. With the typical velocity of the observed MSPs , ~ 70 km/s, it is expected that 
the displacement of the MSPs (or binary system) with the typical age, ^ 100 Myr, becomes 
larger than the size of the Galaxy. With the relatively slow velocity, V ^ 100 km/s, however, 
the MSPs remain bound in the Galactic potential and hence their Galactic distribution does 
not depend on the spin down age. 

The initial rotation period of MSP is related with the history of the accretion process 
after the decaying stage of the magnetic field (Campana et al. 1999; Takata, Cheng and Taam 
2010b). However, the description of the transition from an accretion powered to the rotation 
powered phase is not well understood due to the complexities in the description of the 
interaction between the magnetosphere of a neutron star and its accretion disk (Romanova 
et al. 2009). Furthermore, the true age of the MSP after the supernova explosion in the 
binary system is different from the spin down age of the MSP. These theoretical uncertainties 
make it difficult to obtain the present distributions of the MSP properties from the initial 
distributions. 

In this paper, we assign the "current" pulsar properties for each simulated MSP, instead 
of modeling from the initial distributions; namely, we (1) randomly select the age of the 
simulated MSP up to 10 Gyr, (2) shifts the simulated MSP from its birth location to the 
current location, and (3) assign the parameters of the MSP following the observed distribu- 
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tions. We assign the period time derivative (P) and the stellar magnetic field {Bg) following 
the observed P — Bg distribution (Manchester et al. 2005). From the assigned period time 
derivative and the stellar magnetic field, the current rotation period and the spin down age 
are calculated from 

P_3 = 0.97p2p-2o ms (12) 
and 

r = 1.5 X 10^P_3P2^^ yr (13) 

respectively (Lyne & Graham-Smith, 2006). Here Pg is the stellar magnetic field in units of 
10^ G, P_3 and P_2o s'-re the rotation period in units of 1 millisecond and its time derivative 
in units of 10^^*^, respectively. In fact, above process can provide a consistent Galactic 
distribution of the radio MSPs with the observations (c.f. Takata et al. 2011b). 



2.2 Radio emissions 

Using the empirical relation among the radio luminosity, rotation period, and period time 
derivative, the distribution of the radio luminosity at 400 MHz is expressed by (Narayan & 
Ostriker 1990) 

PL,,, = 0.5X'e\ (14) 

where A = 3.6[logio(L4oo/ < Uqq >) + 1.8] with < L400 >= r/lO^-^^P^/VP^, and L400 
is the luminosity in units of mJy kpc^. Here 77 is a scaling factor to adjust the observed 
distribution, and 77 = 1 (or 0.05) for the CPs (or MSPs). The radio fiux on Earth is given 
by >S'4oo = L^qq/cP, where d is the distance to the pulsar. We scale the simulated 400 MHz 
luminosity to the observational frequency using a typical photon index ~ 2 for the CPs and 
~1.8 for the MSPs, respectively (Kramer et al. 1997, 1998). 

We take into account the beaming of the radio emission. For the CPs, we apply the 
half-angle, which is measured from the magnetic axis, of the radio cone studied by Kijak 
and Gil (1998, 2003), 

ucp ~ r.2Ar]llP~'/\ (15) 
where 

r _ An,,-0.26p0.07 pO.3 
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where P_-i5 is the period time derivative in units of 10~^^, and ugh is the radio frequency 
in units of GHz. For the MSPs, the half-angle does not depend on the frequency and is 
approximately described as (Kramer & Xilouris 2000), 

oJmsp ~ uJoiP/l s)-l/^ (16) 

where Uq is randomly chosen in the range 2.75° — 5.4°. The radio emission can be detected by 
observer with a viewing angle between max(0°,a — Ui) and min(a + Ui, 90°), where i = CP 
or MSP, and a is the inclination angle between the rotation axis and the magnetic axis. 

We use the ten radio surveys (Molongo 2, Green Band 2 and 3, Arecibo 2 and 3, Parkes 
1, 2 and MB, Jordell Bank 2 and Swinburne IL), whose system characteristics are listed 
in table 1 of Takata et al (2011a) and the references therein. To calculate the dispersion 
measure, we apply the Galactic distribution of electrons studied by Cordes & Lazio (2002). 

3 7-RAY EMISSION MODEL 

The pulsar rotation energy, which is thought to be essential energy source of the 7-ray emis- 
sion, can be released by both the current braking torque and the magnetic dipole radiation. 
According to the analysis of the force- free magnetosphere done by Spitkovsky (2006), the 
spin down power depends on the inclination angles as oc 1 + sin^ a; in other words, L^d 
changes only by a factor of 2 with the inclination angle. In the present Monte-Carlo study, 
therefore, we ignore the dependence of the spin down power on the inclination angle, and 
apply the conventional expression that Lgd = 4(27r)^i?^_R^/6c'^P^. 

3.1 Two-layer outer gap model 

The outer gap dynamics is controlled by the pair-creation process, which produces a charge 
distribution in the trans-field direction (Cheng, Ho & Ruderman 1986a,b; Takata, Shibata & 
Hirotani 2004). Wang et al. (2010, 2011) argued that the outer gap should be approximately 
divided into two layers, i.e. the main acceleration region starting from the last-open field 
lines and the screening region lying at the upper part of the gap. In the main acceleration 
region, the charge density is ~10% of the Goldreich- Julian value (Goldreich & Julian 1969), 
and a strong electric field accelerates the electrons and positrons up to the Lorentz factor of 
r ~ 10^'^. The accelerated particles emit several GeV photons via the curvature radiation 
process. In the screening region, the large number of pairs created by the pair-creation 
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process starts to screen out the gap electric field. The curvature radiation from the screening 
pairs produces mainly ~100 MeV photons. 

A simple description of the electric field structure inside the two-layer outer gap is dis- 
cussed in Wang et al. (2010, 2011). We denote x, z and (p as the coordinates along the 
magnetic field line, perpendicular to the magnetic field line in the poloidal plane and in the 
magnetic azimuth, respectively. We expect that the particle number density increases expo- 
nentially near the boundary {z = hi) between the main acceleration and screening regions 
(Cheng et al. 1986a,b), and that the charge density is almost constant in the screening region 
(Hirotani 2006). Hence, we approximately describe the distribution of the charge density in 
the 2;-direction with the step function as follows 

,^ I Pi(a;,0), if ^ 2; ^ /ii(x,0), 
p{r) = i (17) 

\^ P2{x,(f)), if hi{x,(f)) < z ^ h2{x,(l)), 

where, I pi I < \pgj\ < IP2I, z = and z = h2 correspond to the last-open field line and 
the upper boundary of the gap, respectively. For simplicity, we define the boundary be- 
tween main acceleration region and the screening region, i.e. hi, with a magnetic field line 
with such approximation that hi/h2 is constant along the magnetic field line. The present 
model predicts that the charge density in the screening region should be proportional to the 
Goldreich- Julian charge density (Wang et al. 2010). This situation will be satisfied because 
of a lot of the pairs created by the pair-creation process in the screening region. In the 
main acceleration region, because the number density is much smaller than the Goldreich- 
Julian value, its distribution along the magnetic field line does not important for the electric 
field distribution. Therefore, we approximate that p — pcj ~ 9{.z,(t))pGj{f) for both main 
acceleration and screening regions, where 

, \ -9i{^). if O^z^/ii(x,0), 

^(2^,0) = < (18) 

y 5(2(0), if /ii(x,0) < 2 ^ /i2(a;,0). 

We assume that gi> ^ and g2> Q so that \p\ < \pgj\ for the main acceleration region and 
IpI > \Pgj\ for the screening region. 

To obtain the typical strength of the electric field in the gap, we find the solution of the 
Poisson equation for each azimuthal angle (c.f. Wang et al. 2010, 2011), 

— $'(X, 2;,0)|^=fixcd = -47r[p(x, 2,0) - PGj(r)]<)!'=fixcd, (19) 

where $' is the electric potential of the accelerating field. Here we assumed that the derivative 
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of the Potential field in the ^-direction is much larger than that of x-direction, and of 0- 
direction. 

In this paper, we neglect the z-dependence of the Goldreich- Julian charge density, and 
approximate the Goldreich- Julian charge density as Pgj{x, 0) ~ —VIBx/2txcRc (Cheng et al. 
1986a,b), where Vt and Rc are the angular frequency of the pulsar and the curvature radius of 
the field line, respectively. The boundary conditions on the lower {z = 0) and upper {z = /12) 
boundaries are given by 

z = 0, 0) = and z = h2,(j)) = (20) 

respectively. Imposing $' and d^' /dz are continuous at the boundary 2; = /ii, we obtain the 
solution as 

nBxhl{xA)\ -gMz'^ + C,z', for O^2;'^/ii(x,0)//i2(x,0) 

$ (r) = < [21) 

cRc h?2(0)(^"-l) + ^i(^'-l), for /ii(x,0)//i2(x,0) ^ ^' ^ 1 



where 

D.(x,6) 



gihiQii - 2/12) + g2{hi - 
hi 

gihl + g2hl 
hi 



'-2 

and z' = z/h2{x,(j)). The accelerating electric field, £"11 = —d^'/dx, is writted as 

^ nBhl{x,(P) [ -gi{<P)z'' + C,{r)z', for ^ z' ^ /ii(x, 0)//i2(x, 0), 

E\\{r) S (22) 

cRs [ g2mz'''-l) + D^{r){z'-l), for /ii(x, 0)//i2(x, 0) < ^' ^ 1, 

where we used the relations of the dipole field that d{Bh\) / dx ~ 0, dz' jdx = d{z/h2)/dx ~ 
0, d{hi/h2)/dx ~ 0, and approximated that dRc/dx ~ 0. 

On the upper boundary, we anticipate that the total potential field (co-rotational poten- 
tial + non co-rotational potential) in the gap is continuously connected to the co-rotational 

potential field outside the gap. This screening condition is described by 
(9$', 
dz 

This condition gives the relation between {hi, /12) and {gi, 5^2) as 

In this paper, we do not consider the azimuthal distribution of the dimensionless charge 
density gi and g2, because we discuss the general properties of the 7-ray emissions. The 
azimuthal structure will be important for explaining the detailed observed properties, such 
as the existing of the third peak in the pulse profile of the Vela pulsar measured by the 



U=h, = -E^{x,z = h2,(l)) = 0. (23) 
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Fermi, and the energy dependence of the pulse phase of the third peak (c.f. Wang et al. 
2011). 

The typical Lorentz factor of the accelerated particles can be estimated by force balance 
between the electric field and the curvature radiation drag force as 

- ^'^^IeX'\ (25) 



V 2e 

The spectrum of the curvature radiation emitted by the individual particle is written as 
^c(i?„r) = ^^F(x), (26) 
where x = E^/E^ with Ec = ShcT'^ /AirRc and 

POD 

F{x) = X K,/,{Od^, 

where -ft's/s is the modified Bessel function of order 5/3. A 7-ray spectrum measured on 
Earth may be expressed by (e.g. Hirotani 2008) 
dF 1 

^ _ ^ N{r.^P^{E,, mdmA, (27) 

where N ~ \pgj\/g is the particle number density, represents the radius to the emission 
point, from which the emission is measured by the observer, AAi is the area of the calculation 
grid in perpendicular to the magnetic field lines. The integrated energy flux between 100 MeV 
and 300 GeV can be calculated from 

r300GeV 

-^7,100 = / —=^dE^. (28) 

JlOOMeV dE^ 

In the gap, the curvature photons are emitted in the direction of the particle motion, 
which may be described as (Takata, Chang & Cheng 2007) 

v = VpB/B + rxCL, (29) 

where the first term represents the motion along the magnetic field line, Vp is calculated 
from the condition that \v\ = c, and the second term is the co-rotation motion. The emission 
direction measured from the rotation axis (i.e. viewing angle) C, and the pulse phase if) are 
calculated from (Yadigaroglu 1997) 

cosC = «.& 

^ = -i^n - r/Ric ■ n, 

where n = v/v, is the unit vector in the direction of the rotation axis, and ipn is the 
azimuthal angle of the emission direction. 
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3.2 Outer gap geometry 

In this paper, we adopt a rotating vacuum dipole field as the magnetosphere, and we assume 
that a strong emission region extends between the null charge surface (f2 ■ 5 = 0) of the 
Goldreich- Julian charge density and the radial distance r = Ric, where Ric = 271 /cP is the 
light cylinder radius. The electrodynamic studies have pointed out that the gap current 
can shift the inner boundary toward the stellar surface (e.g. Takata et al. 2004). However, 
because it is expected that the curvature radiation below the null charge surface appears 
with a emissivity much smaller than that above the null charge surface (e.g. Hirotani 2006), 
we ignore its contribution to the calculation. 

We define the fractional gap thickness measured on the stellar surface as, 

. M^^. (31) 

where Vp is the polar cap radius. Note because the electric field E\\ is proportional to Bh"^, 
it can be found that E\\ oc 

Zhang & Cheng (1997, 2003) have argued a self-consistent outer gap model controlled 
by the photon-photon pair-creation process between the curvature photons and the X-rays 
from the stellar surface. They estimated the gap fraction as 

CP = ^^^^ - = 5.5(P/ls)^^/^^(5./10^^G)"^/^ (32) 

for the CPs and 

fzc, MSP = 7.0 X 10-2(P/1 msY'/'\BjlO'Gr'/'6rl/\ (33) 

for the MSPs. Here 6rr, is the distance (in units of 10^ cm) from the stellar surface to the 
position where the local magnetic field is comparable to the dipole field, and it will be 
Sr^ ~ 1 — 10 cm. 

We note that Zhang & Cheng (1997, 2003) estimated the gap fraction by a completely 
vacuum electric field E\\ = QBf'^Rf^/cRc. With the same gap fraction, the solution described 
by equation ( 122|) gives an electric field at least a factor of four smaller than that used in 
Zhang & Cheng (1997, 2003). This difference can be important for the typical energy of the 
curvature radiation (Ec), because Ec oc E^^^^. In other words, if we derive the gap fraction 
from the pair-creation condition that ExEc = (mgC^)^, where Ex is the X-ray photon energy, 
the present model predicts a fractional gap thickness larger than that of Zhang & Cheng 
(1997, 2003). In this paper, reducing by a factor of four for the electric field in the model of 
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Zhang & Cheng (1997, 2003), we apply the gap fraction increased by a factor of 4^/'^ ~ 1.8 
from those in equations f p2|) and f p3|l . 

Takata et ah (2010a) proposed that the outer gap model can be controlled by the mag- 
netic pair-creation process taken place near the stellar surface. They argued that half of 
particles created in the gap or injected particles into the outer gap at the outer boundary 
return to the stellar surface, and these returning particles emit ~ 100 MeV photons near 
the stellar surface. A good fraction of 100 MeV photons can make pairs by the magnetic 
pair creation process. They argued that if the magnetic field lines near the surface, instead 
of nearly perpendicular to the surface, are bending side-wards due to the strong local field, 
the pairs created in these local magnetic field lines can have a pitch angle, which is defined 
by the angle between the directions of the particle's motion and of the field line, smaller 
than 90°, which results in an outgoing flow of pairs, and hence the pairs control the size 
of the outer gap. We also note that the outgoing flow may be produced by the Compton 
scattering. With this model, the fractional gap thickness is estimated as 



fm - = 0.8Kpl/^ (34) 

Up 



where K ~ B^i2Sj is the parameter characterizing the local parameters, i.e., -Bm,i2 and S7, 
which are the local magnetic field in units of lO^^G and the local curvature radius in units 
of lO^cm, respectively. By fitting the emission characteristic of the 7-ray pulsars observed 
by the Fermi, they estimated as ~ 2 for the CPs and K ~ 15 for the MSPs. When the 
fractional gap thickness fm is smaller (or larger) than /^c, the magnetic pair-creation (or 
photon-photon pair-creation) process controls the gap thickness. 

We note that the outer gap should only exist between the last-open field lines and the 
critical magnetic field lines that have the null charge points at the light cylinder. Figure [U 
shows the polar angle, which is measured from the magnetic axis, of the polar cap rim (6*^, 
solid line) and of the foot points of the critical field lines (6'c, dashed line) as a function 
of the magnetic azimuth. The maximum gap fraction is defined by fmax = {dp — dc)/dpi 
and is represented by the dotted-line in Figure [H The azimuthal angle = in Figure [H 
corresponds to the plain spanned by the rotation axis and the magnetic axis. In this paper, 
we anticipate that the azimuthal expansion of the active gap is limited as / ^ fmax{4') and 
-90° ^(b^ 90°. 
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Figure 1. Latitudes of the polar cap rim dp (solid line) and of the critical magnetic field lines 6c (dashed line). The dotted line 
shows the maximum gap fraction fmax = {9p — Sc)/Op. The magnetic azimuth = corresponds to the magnetic meridian, 
which includes the rotation and magnetic axises. The results are for the inclination angle a = 45° and the rotation period 
P = 0.1 s. 
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Figure 2. Phase plot of the photons having energies larger than 100 MeV. The results are for a = 60°, B = 3 X 10^^ Gauss 
and fzc = 0.1. 
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4 RESULTS 

In this paper, we present the results of the two-layer outer gap model by using the ratio 
of hi/h2 = 0.95 and the dimensionless charge density in the main region of 1 — (71 = 0.3. 
The value hi/h2 = 0.95 is chosen because Wang et al. (2010) fitted the phase-averaged 
spectra of mature pulsars by using hi/h2 ~ 0.95. The value 1 — gi = 0.3 is slightly larger 
than 1 — gi 0.1 used in Wang et al. (2010), but reproduces a more consistent simulated 
population with the Fermi observations. The gap height /12 and the charge density in 
screening region g2 are calculated from equations ( l3Ti) and ( l24l) . respectively. 

4.1 Dependence on the inclination and viewing angles 

The dependence of the characteristics of the calculated 7-ray spectra on the inclination angle 
and the viewing angle are summarized in Figures EMU Figure [2] shows the phase plot of the 
photons having the energies larger than 100 MeV. The result is for a = 60°, = 3 x 10^^ G 
and /^c = 0.1. 

4-1.1 'y -ray flux 

The left panels in Figures [3] and H] show the 7-ray flux (^ 100 MeV) as a function of the 
viewing angle ( and of the inclination angle a. The vertical line (y-axis) represents the 
fractional 7-ray flux, which is deflned by the flux measured by f^Lsd/d"^ (c.f. Takata et al. 
2011). Here we calculated the rotation period from Bg = 3 x 10^^ G and the gap fraction 
/ = fzc- We can see that the calculated flux tends to decrease as the line of sight approaches 
to the rotation axis, where ( = 0°. In the sky map of Figure [21 we see that the large viewing 
angle {( — )■ 90°) can encounter more intense emission region, whereas the small viewing 
angle ^ 90°) will encounter the less intense region or even miss the emission region. 
With the outer gap geometry running from the null charge surface to the light cylinder, 
the observer with a smaller viewing angle {( <^ 90°) may miss the emissions from higher 
latitudes {z ^ h2/2) and measure the only emission from the lower part of the gap (2; ~ 0). 
Because the accelerating electric fleld vanishes at the lower boundary of the gap {z = 0), 
the Lorentz factor of the accelerated particles and the resultant emissivity of the curvature 
radiation around lower boundary of the gap are signiflcantly decreased. Consequently, a 
smaller viewing angle tends to measure a smaller flux of the curvature radiation. In the left 
panels of Figures [3] and HI we also see the tendency that the decrease of the fractional flux 
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with the decrease of the viewing angle is more gradual for larger inclination angle. This 
is because the emission from the outer gap with larger inclination angle covers more wide 
region of the sky. These dependences of the 7-ray flux on the viewing geometry predict that 
the Fermi has preferentially detected the pulsars with larger inclination angles and larger 
viewing angles near 90° (c.f. section l4.2.2p . 

The right panels in Figure [3] and H] summarize the dependence of the fractional 7-ray 
flux on the gap fractional thickness. We can see that the factional flux decreases as the 
gap fraction increases. This dependence is related with the expansion of the outer gap in 
the azimuthal direction. In the present calculation, we have assumed that the azimuthal 
expansion of the active gap is limited as / ^ fmax{4>) and —90° ^ ^ 90°. As the dotted 
line in Figure [T] shows, fmax acquires a maximum value near the magnetic meridian, where 
= 0°, and tends to decrease as the azimuthal angle deviates from = 0. This implies 
that the width of the active gap in the azimuthal direction narrows with the increase of the 
gap fractional thickness, /. Consequently, the fractional flux tends to decrease as the gap 
fraction increases. 



4.1.2 Cut- ojf energy 

The cut-off energy is calculated as the peak energy in the spectral energy distribution of the 
emission from the main acceleration region. Figures [5] and [6] show the dependence of the cut- 
off energies for the CPs and for the MSPs, respectively, on the viewing geometry. The vertical 
axis represents the cut-off energy measured by ShcTQ/AnRic, where Fq = (S-R^^^-En o/2e)^/^ 
with = pB{Ri^. Figures [5] and [6] show that the cut-off energy decreases with the 
decrease of the viewing angle. As we argued in section I4.1.H a large viewing angle C, ~ 90° 
can encounter the emission from the strong accelerating electric field region at the middle of 
the gap [z ~ /i2/2), whereas the observer with a small viewing angle (C <^ 90°) measures the 
emission from the small electric field region near the lower boundary (2; ~ 0). As a result, 
the cut-off energy in the spectrum decreases with the decrease of the viewing angle. We can 
also find in Figures [5] and [6] that the fractional cut-off energy does not depend much on the 
fractional gap thickness, /. 
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Viewing angle ^ (degrees) 

Figure 3. The 7-ray flux (> 100 MeV) measured on the as a function of the viewing angle. The results are for the canonical 
pulsar with = 3 X 10^^ G and / = fzc- Left:The 7-ray flux for the inclination angle is a = 20° (solid line), 40° (dashed line) 
and 80° (dotted line). The results are for the gap fraction of / = 0.25. RightiThe 7-ray flux for / = 0.1 (solid line), 0.2 (dashed 
line), 0.4 (dotted line) and 0.55 (dashed-dotted line). The results arc for the inclination angle of a = 50°. 

4-1.3 Photon index 

We applied the minimized-x^ method to fit the spectrum between 100 MeV and the cut-off 
energy with a single power low form. Figures [7] and [8] show the dependence of the photon 
index for the CPs and MSPs, respectively, on the viewing geometry. Our model predicts that 
the spectral shape is relatively soft with a photon index p ~ 1.8 — 2 for larger viewing angle 
(C — 7- 90°) and hard with p ~ 1.2 — 1.3 for smaller viewing angle (C ^ 90°). In Figures [7] 
and [HI the transition from soft to hard spectra occurs in a narrow range of the viewing 
angle. In the present two-layer model, the screening region and the main acceleration region 
produce the 7-ray photons with a typical energy of ~ 100 MeV and ~ 1 GeV, respectively. 
For the viewing angle closer to C ~ 90°, because the observed 7-ray radiation consists 
of the emissions from both main acceleration and screening regions, the emerging spectrum 
becomes soft with a photon index of p ~ 1.8 — 2 above 100 MeV. For a smaller viewing angle, 
on the other hand, because the emission of the screening region is missing, the emission from 
the only main acceleration region contributes to the spectrum. In such a case, the emerging 
spectrum has a photon index p ~ 1.2 — 1.3 above 100 MeV, which closes to a mono-energetic 
curvature spectrum. 

4.2 Results of the Monte-Carlo simulation 

In this section, we present the results of the Monte-Carlo simulation, in which we assume 
the birth rates of the CPs and of the MSPs are ~0.015 per year and ~ 9 x 10~^ per year, 
respectively. We also assume that the viewing angle and the inclination angle are randomly 
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Figure 4. The same with Figure [S] but for the miUisecond pulsars with Bs = 3 X 10^ G. 
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Figure 5. Dependence of spectral cut-off energy on the viewing geometry. The vertical line represent the cut-off energy in 
units of (3/47r)/icrg/iJ;c, where Tq = (3i?f^E||_o/2e)l/'* with E||^o = f^^iRlc)- The results are for the canonical pulsar with 
Bg = 3 X 10^^ G and / = f-c- The lines correspond to same cases as Figure|3] 
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Figure 6. The same with Figure[5] but for the millisecond pulsars with B^ = 3 X 10^ G. 
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Figure 7. Dependence of photon index on the viewing geometry. The results are for the canonical pulsar with Bj, = 3 X 10 G. 
The lines correspond to same cases as Figure [3] 
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Figure 8. The same with Figure [3 but for the millisecond pulsars with i?s = 3 X 10^ G and / = f^c- 
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Table 1. Population of simulated radio-selected {Nr) and 7-ray-selected {Ng) CPs for six- month, five-year and ten- year Fermi 
observations. The first line; the results for ten radio surveys listed in table 1 of Takata et al. (2010a). The second and third 
lines are the results with the sensitivities increased by a factor of two and ten, respectively, and the bottom is the populations 
associated with only beaming effects of the radio emission. 
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Table 2. The same with Figure 1, but for the MSPs 
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distributed. For the inclination angle close to a ~ 90°, in which the null charge points are 
located on or very close to the stellar surface, it may cause a numerical error when we search 
the last-open field lines of the rotating vacuum field. To avoid it, we limit the inclination angle 
below a ^ 85°. Because a comparison between the simulated and the observed distributions 
for the various properties (e.g. rotation period and magnetic field) of the 7-ray pulsars was 
done in Takata et al. (2011a,b), we focus on the characteristics of the 7-ray radiations with 
the viewing geometry. 

For the sensitivity of the observations, we refer the Fermi first pulsar catalog (Abdo et al. 
2010a). For the radio-selected pulsars, the Fermi archived the sensitivity F ~ 10~^^ erg/cm^s 
(or ~ 3x 10"^^ erg/cm^s) for the galactic latitudes \b\ ^ 5° (or < 5°) with the six-month long 
observations. For the CPs, the Fermi six-month data allow us to detect the pulsed period 
by blind search, if the 7-ray fiux is larger than F ^ 2 x 10~^^ for \b\ ^ 5° and F ^ Q x 10"^^ 
for |6| ^ 5°. For the MSPs, because there are no such detections of the 7-ray-selected pulsar 
so far, we cannot simulate the Fermi sensitivity of the blind search. In this paper, therefore, 
we simulate the population of the 7-ray-selected MSPs with the Fermi sensitivity of the 
blind search of CPs. 

4-2.1 Viewing geometry 

Figures [9] and [TO] plot the inclination angle (a) and the viewing angle {() for the simulated 
100 canonical and millisecond 7-ray pulsars, respectively. The fiUed-circles and the boxes 
represent the radio-selected and 7-ray-selected 7-ray pulsars, respectively. For the CPs, we 
can see in Figure M that the radio-selected and 7-ray-selected pulsars distribute at different 
region in a — C plane. Specifically, the radio-selected CPs group together around the line 
a = (. This is because we have assumed that the magnetic axis is the centre of the radio 
cone, indicating that the radio emission can be detected by observer with a viewing angle 
( a. The scattering from the line a = C is related with the width of the radio cone. 
For example, if we assume a narrower cone than that of equation f|T5|) . the amplitude of 
the scattering is reduced. Figure [9] also shows that the detected 7-ray pulsars are mainly 
distributed with a viewing angle C ~ 70 — 90°. These results are consistent with the results 
of Watters & Romani (2011). 

Unlike with the CPs, the distributions of the radio-selected and 7-ray-selected MSPs 
on the a-( plane are overlapped each other. This is because the width of the radio beam 
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described by equation (fT6l) covers almost whole sky. We note that most 7-ray-selected MSPs 
in the simulation irradiate the Earth with the radio emissions, but the radio fluxes are lower 
than the sensitivity of the simulated radio surveys (c.f. Tables 1 and 2). 

In Figures |9] and [THl we can see that no 7-ray pulsars are detected with a smaller incli- 
nation angles {a ^ 90°) and a smaller viewing angles (C ^ 90°). This is because the 7-ray 
flux decreases as the viewing angle and/or the inclination angle decrease, as we discussed in 
section l4rT] (c.f. Figures [3] and S]) . Hence, our simulation results predict that the pulsars with 
larger inclination and larger viewing angles (a, — )■ 90°) have been preferentially detected 
by the Fermi six-month long observations. 

4.2.2 Population 

Tables 1 and 2 show the simulated population of the CPs and MSPs, respectively, with the 
Fermi six-month, five-year and ten-year long observations. Here we scale the sensitivity of 
the Fermi observations as oc -\/T, where T is the length of the observation time. In addition, 
the second lines ("Ra. Sen. (x2)") and the third lines ("Ra.Sen. (xlO)") in the tables show 
the results for the ten radio-surveys but we increase the sensitivities by a factor of two and 
of ten, respectively, and the fourth lines ( "Beaming" ) show the population associated with 
the only beaming effect of the radio emission. 

With the previous radio surveys (first line in Table 1), the present simulation shows 
that 40 radio-selected and 39 7-ray-selected CPs can be detected by the Fermi six-month 
observations, indicating that the present model predicts more 7-ray pulsars than the Fermi 
observations (~ 20 for both radio-selected and 7-ray-selected 7-ray pulsars). To explain the 
difference between the simulated and observed numbers, the several reasons will be expected; 
(1) the Fermi sensitivity will become much worse at the Galactic plane, and (2) the 7-ray 
emissions from the pulsars will be missed by source confusion with the complex regions 
and the unresolved sources that are not modeled in the diffuse backgrounds. The predicted 
number 10 of the radio-selected 7-ray MSPs is consistent with the observed number 9. The 
present model predicts 56 (or 59) radio-selected and 138 (or 182) 7-ray-selected CPs can 
be detected by the Fermi with five-year (or ten-year) observations. For the MSPs, 14 (or 
16) radio-selected 7-ray pulsars will be detected by the Fermi with five-year (or ten-year) 
observations. 

We see in the first lines of Tables 1 and 2 that the simulated numbers of radio-selected 
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7-ray CPs and MSPs increase only ~20 and ~10 sources, respectively, over even ten-year 
Fermi observations. This implies that most presently known radio pulsars (~ 2000 for CPs 
and ~ 80 for MSPs) might not be discovered by the Fermi. For the 7-ray-selected pulsars, 
on the other hand, the simulation predicts that the Fermi can detect about 140 CPs for 
about five-year observations. We note that the predicted number of the radio-loud (or radio- 
quiet) 7-ray pulsars really depends on the sensitivities of radio surveys, as Tables 1 shows; 
for example, the predicted number of the 7-ray-selected CPs after ten-year of the Fermi 
observations decreases from 182 to 120 if the sensitivities of the radio surveys increases by a 
factor of ten. This indicates that a deep radio search may find more radio-emissions from the 
7-ray-selected pulsars, such as LAT PSRs J1741-2054 and J2032+4127 (Camilo et al. 2009). 
Table 1 also shows if the sensitivities of the radio surveys increase by a factor of ten, the 
ratio of the radio-loud and radio-quiet 7-ray pulsars is almost determined by the beaming 
effect of the radio emissions. 

For the MSPs (Table 2), most simulated pulsars are categorized as the 7-ray-selected 
pulsars with the previous sensitivities of the radio surveys, although the Fermi has not 
confirmed the radio-quiet MSPs. We argue that it may be difficult to identify radio-quiet 
MSPs, because the detection of the rotation period by the Fermi blind search is much 
harder than that of the CP. Furthermore, if the MSP is in binary system, the effects of 
the orbital motion on the observed rotation period make it even harder to confirm the 
millisecond rotation period by the blind search. Takata et al. (2011b) discussed that the 
7-ray-selected MSPs in the simulation correspond to the Fermi unidentified sources located 
at higher Galactic latitudes. 

We note that the radio cones from MSPs are quite huge so that most MSPs irradiate 
the Earth with the radio emissions, as the bottom line in Table 2 shows. This implies that 
more radio MSPs associated with the Fermi unidentified sources will be detected by the 
future radio surveys, such as the discovery of new 20 radio MSPs associated with the Fermi 
unidentified sources (Ray 2010; Caraveo 2010; Ransom et al. 2011; Keith et al. 2011). As 
third line in Table 2 shows, however, even we drastically increase the radio sensitivity by a 
factor of 10, the number of radio-selected MSPs detected by 10 year Fermi observations can 
increase from 16 to 94, it is still much less than the expected 227 7-ray-selected millisecond 
pulsars. Unless the Fermi sensitivity of the blind search is improved, the most 7-ray MSPs 
will not be identified and will contribute to the Fermi unidentified sources and/or the 7-ray 
background radiations. 
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4-2.3 ^-ray radiation characteristics 

Figures [TTIfT3] show the radiation characteristics (7-ray luminosity, cut-off energy and photon 
index, respectively) versus the characteristics (the spin down power or the magnetic field 
strength at the light cylinder) for the 7-ray pulsars with the six-month observations. In 
those figures, we plot the Fermi data with errors taken from Abdo et al. (2010a) and 
Saz Parkinson et al. (2010), and present the simulated pulsars in the sub-figures. For the 
simulated pulsars, we randomly choose 200 simulated 7-ray pulsars (except for the 7-ray- 
selected MSPs). 

In Figure [HI the 7-ray luminosity is calculated from = Aird^Fj^ioo, where d is the 
distance, and the histograms of the values of the 7-ray luminosity are projected along the 
right-hand axis. The solid and dashed histograms represent the distributions for the simu- 
lated and observed 7-ray pulsars, respectively. The present model predicts that most 7-ray 
CPs have a spin down power of Lsd ~ io^^~'^^ erg/s and a 7-ray luminosity of ~ 10^^"'^^, 
while MSPs have a Lgd ~ IQ^'^"^^ erg/s and ~ ^g32. 5-34.5 erg/s, which are consistent with 
the Fermi observations. 

In Figure [HI the spin down power Lgd and the 7-ray luminosity L-y of the simulated 
pulsars can be related as oc Lf^ with /3 ~ for Lsd ^ iQ^s-se erg/s and f3 ~ 0.5 for 
Lsd ^ 10^^"'^^ erg/s. In the present emission model, the 7-ray luminosity is proportional to 

oc f^Lsd as Figure [3] and [H indicate (c.f. Takata et al. 2011b). The change of the slope 
is caused by switching gap closure process between the photon-photon pair-creation process 
and the magnetic pair-creation process. As the equations (j32l) -( jMl) show, the gap fraction 
depends on the rotation period and the magnetic field as fzc oc p26/2i^-4/7 ^j^g photon- 
photon pair-creation process and oc P^/^ for the magnetic pair-creation process. These 
relations imply that the 7-ray luminosity depends on the spin down power as oc L^.^J'^ 

5 /8 

for the photon-photon pair-creation process and oc LJ^ for the magnetic pair-creation 
process. Equating fzc ( j32l) and fm (j34l) corresponds to Lsd ~ 10^^"'^^ erg/s. The change of 
the slope /3 has been found by Wang et al (2010), who used the two-layer outer gap model 
to fit the phase-averaged spectrum of the mature 7-ray pulsars observed by the Fermi. 

Although the general trend of the relation between the 7-ray luminosity and the spin 
down power is explained with a simple form oc Lf^, we can see in Figure [H] that some 
simulation samples deviate from the relation. For example, some simulated pulsars with a 
spin down power of Lsd ~ 10"^^ erg/s have a 7-ray luminosity of L^ ~ ]^g32-33 erg/s, which 
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is about two or three order smaller than the typical value of ~ 10^^ erg/s. We emphasize 
that this low efficiency of the 7-ray emission is mainly caused by the effects of viewing angle. 
The pulsars lying on the relation oc L^^ are observed with viewing angles of C ~ 90°, 
whereas those ("apparently") low-efficient 7-ray pulsars are observed with smaller viewing 
angles. Although the flux depends on also the inclination angle, the flux is more sensitive to 
the viewing angle than the inclination angle, as Figures [3] and H] show. We note that lower 
efficient 7-ray pulsars tend to locate closer to the Earth. 

Figure [12] shows the cut-off energy versus the magnetic field at the light cylinder, Bi^. 
The symbols and the histograms correspond to same cases as Figure [Til As sub-figure in 
Figure [I2] indicates, the present simulation predicts that most 7-ray CPs and MSPs have a 
magnetic field at the light cylinder of Bi^ ^ lO"' G and Bi^ ^ 10"^ G, respectively, and have 
a cut-off energy smaller than ~ 2 GeV. These features will be consistent with the Fermi 
observations. In the figure, some simulated CPs have a cut-off energy significantly smaller 
than the typical value E'c ~ 2 GeV. As well as the case of the 7-ray luminosity, this deviation 
from the typical value is cause by the effects of the viewing angle, as Figures [S] and [HI imply. 
In Figure [121 our simulation predicts that the typical cut-off energy (1-1.5 GeV) of the MSPs 
is smaller than that (~ 2 GeV) of the CPs. It is difficult to discuss the difference in the 
observed cut-off energies between the CPs and the MSPs, because of the large observational 
errors. 

Figure [T3l represents the photon index versus the spin down power. We can find that the 
model distribution of the photon index has two peaks at p ~ 1.2 — 1.3 and p ~ 1.8 — 2, 
and the observed distribution (dashed histogram) may also have two peaks at p ~ 1.3 and 
p ~ 1.7. With the present model, the hard component ofj»~1.2 — 1.3 corresponds to the 
spectrum associated with the emission from the only main acceleration region, while the soft 
component p ~ 1.8 — 2 corresponds to the spectrum composed of the emissions from both 
main and screening regions, as we discussed in section 14.1.31 The present model does not 
predict very hard spectrum with a index p ^ 1, which has been indicated for some Fermi 
pulsars. 

Figures [H] and [161 represent the averaged apparent fractional thickness, which is defined 
by fa = {AncfF^/LsdY^^, and the photon index, respectively, as a function of the spin down 
power, and Figure [T5l shows the averaged cut-off energy as a function of the magnetic field at 
the light cylinder. In the figures, the solid and dashed lines represent the results for the radio- 
selected and 7-ray-selected CPs, respectively. In addition, the thick and thin lines correspond 
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to the results of the simulation and the Fermi observations, respectively. In Figure [TU our 
model predicts a tendency that the apparent fractional thickness, /a, tends to decrease with 
the increase of the spin down power. This is because the true fractional thickness, fzc or 
fmi tends to decrease with the increase of the spin down power, that is, f^c L'J^^'^^ B^^'^^ 
from equation (|32|) and fm oc L^^^^B^^'^ from equation flM|l . Because the photon index of 
the spectrum tends to increase with the decrease of the fractional gap thickness as Figure [7] 
shows, the averaged photon index in Figure [T6] increases with the spin down power. We find 
that these behaviors are qualitatively consistent with the Fermi observations. In Figure [T5l 
our model predicts that the typical cut-off energy does not depend much on the magnetic 
field strength at the light cylinder, while the present Fermi data may have a tendency that 
the cut-off energy increases with the magnetic field strength at the light cylinder. However, 
because the observational errors are so large, a more deep observation to reduce the errors 
may be required to discuss the tendency. 

In Figures [171 and [18], we summarize the distributions of the radiation characteristics for 
the simulated canonical and millisecond 7-ray pulsars, respectively, including both radio- 
selected and 7-ray-selected pulsars. The solid and dashed lines are results for the simulated 
six-month and ten-year Fermi observations, respectively. Comparing the solid and dashed 
lines, we find that the distributions do not depend much on the time span of the Fermi 
observations. However, we note that a longer observation enables to detect 7-ray pulsars 
with smaller fluxes, which include 7-ray pulsars with the viewing angle close to the rotation 
axis (c.f. Figures [3] and H]) . For the observer with a viewing angle close to the rotation axis, 
the photon index tends to be p ~ 1.2 — 1.3, as Figures [7] and [H] show. Therefore, our model 
prediction is that a longer Fermi observation will detect more 7-ray pulsars with photon 
indexes p ~ 0.12 — 0.14, as right panels in Figures [T7I and [T8] indicate. 

4.2.4 Pulse profiles 

Figures [T91I22] present the calculated pulse profiles for the different type of the 7-ray pulsars. 
For each type of the 7-ray pulsars, 64 samples are randomly chosen to present the model 
prediction on the statistical distribution of the morphology of the pulse profiles. From the 
left to right panels and from the top to bottom panels in the figures, the inclination angle 
increases. In principle, we can quantitatively compare the simulated morphology of the 
pulse profiles in Figures [T91I22] (e.g. phase-separation, number of peak) with the Fermi 
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observations. However, we would like to point out that it is very difficult to quantify the 
morphology to compare with the observations. In Figure [191 for example, the pulse profile 
represented in upper-left conner is indeed a double peak structure with the narrow phase- 
separation between two peaks. In the observations, however, the identification of the double 
peak structure with such narrow peak separation really depends on source counts and timing 
ephemerides of the pulsars. In Figure [20| the simulated pulse profile presented at the seventh- 
line and the ffist-column [a = 74.2° and ( = 67.5°) has the ffist peak much smaller than 
the second peak. In the observations, the detection of such a small peak will depend on 
the strength of the background radiation. Moreover, because the intensity of pulse peak will 
depend on the modeling of the gap structure, it is difficult to discuss the distribution of the 
intensity ratio of the first and the second peaks with the present simple three-dimensional 
model. In the present paper, therefore, we avoid a quantitative discussion for the morphology 
of the pulsed profile. Instead, we can only present a qualitative comparison. 

We note that the peaks emerging in the calculated pulse proffies are caused by the so- 
called caustic effect (Romani and Yadigaloglu 1995; Cheng, Ruderman and Zhang 2000; 
Dyks, Harding and Rudak 2004), in which more photons are observed at narrow width of 
the rotation phase due to the special-relativistic effects, i.e. the aberration of the emission 
direction and photon's travel time. In such a case, the pulse phase and the number of the 
main peak are not affected much by modeling of the gap structure, whereas the intensity of 
the peak depends on the gap structure. 

We can see in Figures fT9] and [20] that the pulse profiles of the simulated CPs (in particular, 
7-ray-selected CPs) are described by double peak structure rather than single peak profile. 
This result may explain the tendency that most canonical 7-ray pulsars observed by the 
Fermi show the double peak structure. In the present simulations, we have predicted that 
the 7-ray pulsars measured from ( ~ 90° are preferentially detected in the simulation (c.f. 
section [¥?T]) . With the viewing angle closer to C ~ 90°, the calculated pulse proffie tends to 
have the double peak structure. 

For the radio-selected canonical 7-ray pulsars, although the double peak structure are 
more or less common feature, single pulse structure or the double peak structure with a 
narrow phase separation stands out for pulsars with inclination angles a ~ 40° — 50°, as we 
can see in Figure [191 This is because the simulated radio-selected canonical pulsars have the 
viewing angle similar to the inclination angle (c.f. Figure [9}), ( ^ a. With the present outer 
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gap geometry, the pulse profile of the smaller viewing angle {C, ^ 90°) tends to have single 
peak or double-peak with a narrow phase separation, as the phase plot of Figure [2] shows. 

Comparing the pulse profiles of the CPs (Figures [T9l and [201) and of the MSPs (Figures [21] 
and [22]) . one may see that the pulse profile with single peak is more common for the MSPs 
than the CPs. In the present model, the fractional gap thickness fzc,MSP of equation ( 133|) 
for the MSPs tends to be larger than fzc,cp of equation ([321) for the CPs. With a larger 
fractional gap thickness, the emission from the higher altitude (larger z in units of the 
light radius) contributes to the pulse profile. Because the emission from the higher altitude 
produces two caustic peaks with a narrower phase separation, the calculated pulse profiles 
of the MSPs have single peak or double peak structure with the narrow phase separation 
more than that of the CPs. 

5 SUMMARY AND DISCUSSION 

In this paper, we have applied the so called two-layer outer gap model for the 7-ray radiations 
from the pulsars, and have focused on the dependence of the 7-ray radiation characteristics 
on the inclination angle and the viewing angle. We showed that the 7-ray flux and the 
spectral cut-off energy decreases as the viewing angle deviates from ( = 90°. The spectrum 
above 100 MeV becomes soft with a photon index p ~ 1.8 — 2 for the observer with a 
larger viewing angle {( — 90°), whereas it becomes hard with a index p ~ 1.2 — 1.3 for 
the observer with a smaller viewing angle {( ^ 90°). The spectrum with a photon index 
p ~ 1.8 — 2 consists of the emissions from both main acceleration and screening regions, 
while p ~ 1.2 — 1.3 corresponds to the emission from the only main acceleration region. 

We have developed the Mote-Carlo simulation for the population of 7-ray emitting pul- 
sars. The our simulation predicts that 56 (or 59) radio-selected, 138 (or 182) 7-ray-selected 
canonical 7-ray pulsars and 14 (or 16) radio-selected 7-ray MSPs can be detected by five- 
year (or ten-year) Fermi observations. Even we drastically increase the radio sensitivity by a 
factor of ten, the most simulated MSPs are expected as the 7-ray selected MSPs, which will 
contribute as the Fermi unidentified sources and/or 7-ray background radiations. For the 
viewing geometry, our simulation predicts that the radio-selected CPs have the inclination 
angle (a) and the viewing angle {() of a ~ ^^'^^ that most 7-ray pulsars detected by the 
Fermi have been measured with a viewing angle of C = 70° ~ 90°. We demonstrated that 
the spin down power Lg^ and the 7-ray luminosity of the simulated pulsars is related 



Population study for ^-ray pulsars 29 

as oc Lf^ with /3 ~ for L^d ^ lO^^-se g^g/s and /3 ~ 0.5 for L^d ^ 10^^~^^ erg/s. The 
change of the slope (3 is associated with the switching of the gap closure mechanism from the 
photon-photon pair-creation process of Lgd ^ igss-se g^g^'g ^j^g magnetic pair-creation 
process of Lgd ^ io^^~^^ erg/s. We showed that the distribution of the photon index of the 
simulated 7-ray pulsars has two peak at p ~ 1.8 — 2 and p ~ 1.2 — 1.3, which may explain 
the observed distribution. We expect that more 7-ray pulsars with the hard spectrum of 
p ~ 1.2 — 1.3 will be detected by the future Fermi observations. For the pulse profiles, the 
present simulation will explain the observational tendency that most canonical 7-ray pulsars 
detected by the Fermi have the pulse profiles with the double peaks. The our model expects 
that single pulse profile or the double pulse profile with a narrower phase separation is more 
common for MSPs than CPs. 

The present model can be used to diagnose the viewing geometry of the 7-ray pulsars. 
For example, the present model predicts that most present 7-ray pulsars detected by the 
Fermi will have the viewing angle, ( ~ 70° — 90°, and hence the pulse profile with the 
double peaks is more common than that with the single peak. On the other hand, we expect 
that PSR J0659-I-1414, which is known as the extremely low efficient 7-ray pulsar (Abdo et 
al. 2010c), has a unique viewing geometry. PSR J0659+1414 is the radio-loud 7-ray pulsar 
and its spin down power is Lsd = 3 x 10'^'^ erg/s. Interestingly, the 7-ray emissions from 
PSR J0659+1414 has been observed with a 7-ray luminosity, ~ 3 x 10^^ erg/s, which is 
about two order smaller than the typical value for the pulsars with a similar spin down power. 
To explain the observed low efficiency of the 7-ray radiation, our simulation predicts that 
the inclination angle and the viewing angle are relatively small, say a ~ C ~ 40° — 50°, and 
hence the apparent efficiency is significantly reduced from the typical value (c.f. Figure [3l). 
PSR J0659+1414 is also known that the spectral cut-off energy, Ec ~ 0.7 GeV, is suggestively 
smaller than the typical value of -Ec ~ 2 GeV. This behavior is also able to be explained 
by the effects of the viewing geometry, as Figure [5] indicates. Furthermore, our simulation 
expects that the pulsars measured as the low efficient 7-ray pulsars are located closer to 
Earth. In fact, PSR J0659+1414 with d ~ 0.28 kpc is one of the nearest 7-ray pulsars. 
Finally, the observed broad single pulse profile of PSR J0659+1414 is also expected by the 
present calculation. As discussed in section 14.2.41 (Figure fT9|) . the single pulse structure or 
double peak structure with a narrow phase separation is more common for the radio-selected 
canonical 7-ray pulsar with an inclination angle a ~ 40° — 50°. On these ground, the unique 
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radiation properties of PSR J0659+1414 can be explained, if the Fermi has measured the 
7-ray emission from the our gap with the viewing geometry a ~ ~ 40° — 50°. It is expected 
that the population of the low-efficient 7-ray pulsars will be increased by the future Fermi 
observations. The diagnose its 7-ray efficiency, pulse profile and cut-off energy will provide 
a strong constraint on the emission models. 
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Figure 9. The rotation axis a and the viewing angle for the simulated canonical 7-ray pulsars. Filled-circles: Radio-selected 
pulsars. Boxes:7-ray-selected pulsars. 
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Figure 10. The same with Figure [S] but for the miUisecond pulsars. 
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Figure 11. The 7-ray luminosity, = i-wd? F^^iqq, versus the spin down luminosity, L^^i, as result of th six-month long 
observation. The Fermi data with errors are taken from Abdo et al. (2010a) and Saz Parkinson et al. (2010), and the 200 
samples of the simulated pulsars are plotted in the sub-figure, in the sub-figures. Filled-circles;Radio-selected CPs. Square:7- 
ray-selected CPs. Filled-trianglesiRadio-selected MSPs. The histograms of the values of the 7-ray luminosity are projected 
along the right-hand axis. The solid and dashed histograms are results for simulated and observed pulsars, respectively. 



Population study for ^-ray pulsars 35 



10 



> 

CD 
O 

>^ 
D5 

05 1 
a> 

I 

o 



0.1 



if 



Radio-selected CPs« 
y-selected CPsn 
Radio-selected MSPsa 



10 



0.1 



i i 




3.5 



4.5 



5.5 



6.5 



logioBic[G] 



Figure 12. The cut-off energy versus the magnetic field at the hght cyhnder, Bic- The symbols and histograms correspond to 
same cases as Figure [TT] 
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Figure 13. The photon index versus the spin down luminosity. The symbols and histograms correspond to same cases as 
Figure [TT] 
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Figure 14. Averaged value of the observed fractional gap thickness, fo = {^ircP F.y / L gj^)^ ^ ^ , of the detected 7-ray pulsars as 
function of the spin down power. The solid and dashed histograms show the results of the radio-selected and 7-ray selected 
pulsars, respectively. The thick and thin lines correspond to the simulation and Fermi observation, respectively. 
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Figure 15. Averaged value of the cut-off energy of the detected 7-ray pulsars as function of the spin down power. The 
histograms correspond to same cases as Figure [HI 
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Figure 16. Averaged value of the photon index of the detected 7-ray pulsars as function of the spin down power. The histograms 
correspond to same cases as Figure [Ml 
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Figure 17. The simulated distributions for tlic canonical pulsars. Loft;Luminosity. Middlc:Cut-off energy. RightiPhoton index. 
The solid and dashed lines arc results for simulated six-month and ten-year Fermi observations, respectively. 
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Figure 18. The same with Figure [TTI but for the millisecond pulsars. 
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Figure 19. The pulse profiles of 64 samples of the simulated radio-selected 7-ray pulsars. 
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Figure 20. The pulse profiles of 64 sample of the simulated 7-ray selected canonieal pulsars. 
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Figure 21. The pulse profiles of 64 sample of the simulated radio-selected millisecond 7-ray pulsars. 
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Figure 22. The pulse profiles of 64 sample of the simulated 7-ray- selected millisecond pulsars. 



